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Abstract 

Double-7 vibrations in deformed nuclei are studied in the context of the 
interacting boson model with special reference to their anharmonic character. 
It is shown that large anharmonicities can be obtained with interactions that 
are (at least) of three-body nature between the bosons. As an example the 7 
vibrations of the nucleus 1 ggEr 98 are studied in detail. 
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Nuclear quadrupole shape oscillations can be of two types: [3 or 7 vibrations J]J . The (5 
vibration preserves axial symmetry and a one-quantum excitation gives rise to a K = band 
where K is the projection of the angular momentum on the axis of symmetry of the nucleus. 
A 7 vibration breaks axial symmetry and leads to a K = 2 band. Although their existence 
has been conjectured a long time ago the observation and interpretation of /3-vibrational 
K n = + bands is still fraught with questions and difficulties. In contrast, 7- vibrational 
K w = 2 + bands are systematically observed in deformed nuclei and their properties are 
correspondingly better understood. 

Since single-7 vibrations are so well established, it is natural to search for double-7 
vibrations and to examine their harmonic nature (i.e., whether they occur at twice the 
energy of the single vibration.) Two intrinsic K = 2 quanta can be combined parallel or 
antiparallel and hence lead to two bands: one with K = and another with K = 4. The 
experimental identification of double-7 vibrations in deformed nuclei is difficult since they are 
expected to lie above the pairing gap and to mix with two-quasiparticle excitations, resulting 
in fragmentation and a corresponding reduction in the collectivity of the states. During the 
last few years, however, a steady improvement of experimental techniques has allowed the 
measurement of low-spin states in the energy region of interest This possibility has 

reopened the old debate on the existence of two-phonon {(3 or 7) vibrational states and 
their properties. Experiments have been reported recently pointing out the existence of 
double-7 vibrations in several deformed nuclei with a wide range of anharmonicities |6|-[TT 



In particular, in Refs. the first observation of the K n = + and K 1 * = 4 + double-7 
states in one nucleus, 166 Er, is reported. They are observed at 1.949 MeV and 2.029 MeV, 
respectively. This information is of great interest since it provides a stringent test of nuclear 
models; for instance, the Quasiphonon Nuclear Model (QPNM) predicts no K v = + two- 
phonon state below 2.5 MeV in 166 Er |I2"|| . Several calculations of two-phonon states, using 
either phenomenological or microscopic models, are available, particularly for 166 Er and 168 Er 
|T2^|TT||. One of the models employed is the Interacting Boson Model (IBM) In the 

simplest version of this model, referred to as IBM-1, an even-even nucleus with n valence 



Tl 

nucleons is treated as a system of N = — bosons with £ = (s bosons) or t = 2 (d bosons). 
In the usual formulation of the model only up to two-body interactions between the bosons 
are taken. 

What are the predictions of IBM with regard to two-phonon states in deformed nuclei and 
their (an)harmonic nature? It was pointed out some time ago by Bohr and Mottelson [13j 
that the IBM-1 is unable to accommodate large anharmonicities, as observed for instance 
in 168 Er. Subsequently, it was shown that these can be described but require g bosons with 
I = 4 in addition to the s and d bosons (sdg-lBM) More recently, we reported a study 
of two-phonon states in IBM-1 treated in the intrinsic frame fll| and showed that the IBM-1 
is a harmonic model in the limit of large boson number. Anharmonicities can only exist 
for finite boson number and they are always small if only up to two-body interactions are 
considered. It was also suggested that anharmonicity in the model is linked to triaxiality. 
Since it is known that IBM-1 with only up to two-body interactions cannot give rise to a 
stable triaxial minimum, the model's capability for describing anharmonicities depends on 
the inclusion in the Hamiltonian of higher-order interactions, some of which are known to 



induce triaxial shapes [20.21]. 



In this article the relation between three-body interactions in IBM-1 and the anhar- 
monicity of 7 vibrations in deformed nuclei is investigated. Although the analysis presented 
is not exhaustive, it is shown that anharmonic behavior can be obtained with reasonable 
three-body interactions. As an example, the energy and E2 transition properties of the 7 
vibrations of the nucleus 166 Er are studied in detail. In addition, the nature of the OjJ" state 
in the same nucleus, which has been the subject of an intense debate in the last few years 
[22]-p5|, is briefly discussed. 



The Hamiltonian adopted in the following includes a quadrupole-quadrupole term, a 
rotational L 2 term, and three-body interactions between the d bosons, 

H = -kQ -Q + kL-L + Y^Oi (V x rft)( fc ) x ■ ((d x t£) (fc) x d) (l) , (1) 



kl 

where • denotes scalar product, d^ = (— l) M oL M , Q is the boson quadrupole operator, and L 



is the angular momentum operator: 

7\ ( 2 ) f /77T / ,+ r\ (!) 



Q = s^d+(tfs + x((tf x dy ', L = VTo (d ] x d) . (2) 

Five independent three-body ci-boson interactions exist which have Z = 0, 2, 3, 4, and 6. 
Interactions with the same I but different k are not independent but differ by a normalization 
factor only p0| . The combinations (k,l) = (2,0), (0,2), (2,3), (2,4), and (4,6) are chosen 



here. 

The Hamiltonian (H) is certainly not the most general that can be considered. Notably, 
a vibrational term e^d which dominates in spherical nuclei is omitted since it is thought of 
lesser importance in the deformed nuclei considered here. It is clear that the inclusion of 
such additional terms might improve the quality of detailed fits to particular nuclei such as 
the one for 166 Er presented below. Finally, of all possible three-body interactions only those 
between the d bosons are retained here since these are most crucial for obtaining a stable 



triaxial minimum ||20|| . 



For the discussion of anharmonicities of 7 vibrations it is convenient to define the fol- 
lowing ratios of excitation energies: 



£x(0+ ) _ ff x (4+ 7 ) - E x (4j) 

° " £ x (2+) - E x (2f) ' 4 " £x(2+) - E x (2f) ' 1 } 

where 0+ 7 and 4+ 7 are the band heads of the K n = + and K n = 4 + double-7 bands, 
respectively. It should be noted that the quantities FO K do not depend upon the L 2 term 
in the Hamiltonian; if a single three-body term is included they depend on two parameters, 
X and the ratio 81/ k. In the present work the identification of the states 77 and 4+ 7 is 
based on the B(E2) values for decaying to the single gamma state. In Fig. [I] the quantities 
FC K are plotted (for iV = 15 bosons) as a function of the quadrupole parameter \ (varying 
between its SU(3) and 0(6) values — |v^7 and 0) in the absence of three-body interactions. 
The ratio Rj remains about constant and of the order 1.8; Rq shoots up for small |x|. Close 
to the 0(6) limit the concept of a 7 vibration is not well defined and so nothing is plotted 
for \x\ < 0.15. The value of x is constrained by E2 transition probabilities and in deformed 



rare-earth nuclei it ranges typically between —0.4 and —0.7 | 26| . From Fig. [I] it is clear that 
no substantial anharmonicity occurs in the 7 vibration for these values of \- 

In Fig. the influence of the various three-body interactions is shown for a typical value 
°f X (% = — 0.5) and for N = 15 bosons. It is seen that 7-vibrational anharmonic behavior 
is obtained which can be different for the K 71 = + and K 71 = 4 + bands (e.g., positive for 
the former while negative for the latter.) Care has been taken to plot results only up to 
values of 61 that do not drastically alter the character of rotational spectrum; beyond these 
values, the three-body interaction, being of highest order in the Hamiltonian (|]), becomes 
dominant. Also shown in Fig. [| are the ratios tV K as observed in 166 Er |7|||, Rq = 2.76 
and Rj = 2.50. This simple analysis shows that, on purely phenomenological grounds, the 
appropriate three-body interaction with the correct anharmonic character for the = + 
and K* = 4+ bands in 166 Er, has I = 4. 

Figure || shows the experimental spectrum of 166 Er [7|§ and compares it to the eigenspec- 
trum of Hamiltonian ([]]) with an I = 4 three-body interaction. The parameters are k — 23.8 
keV, x = — 0.55, k' = —1.9 keV, and #4 = 93.9 keV, with boson number N = 15. With 
these values the calculated excitation energies of the double-7 band heads are 1926 keV and 
1972 keV for the K n = + and K 71 = 4 + levels, respectively, leading to the ratios Rq = 2.82 
and Rj = 2.45, in excellent agreement with observation. Note, however, that although all 
7-band heads are well reproduced by the calculation, problems arise for the moments of 
inertia, in particular of the 7 band. An extensive survey of combinations of cubic <i-boson 
interactions has shown that it is difficult to substantially improve upon this fit although it is 
of course near-impossible to do an exhaustive search of the complex parameter space of all 
three-body interactions. In contrast, exploratory searches with simple quartic Hamiltonians 
quickly yield the correct result with respect to both band-head energies and moments of 
inertia. 

For the calculation of E2 transition probabilities the Consistent-Q Formalism (CQF) |2j 
is adopted by using the E2 transition operator 
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f (E2) = e cS Q, (4) 

where Q is the boson quadrupole operator used in the Hamiltonian (0) and e e fj is a boson 
effective charge, determined from the observed B(E2;2± — > 0+) value. It should be noted 
that the inclusion of three-body terms in the Hamiltonian would allow the use of a two-body 
E2 operator. However, we have not tried to do that in order to keep the calculation on the 
anharmonicity of the double-gamma excitation as simple as possible. In table | the observed 
B(E2) values and ratios concerning the 7 vibrational band heads in 166 Er are summarized 
and compared to the theoretical results obtained with e^ s = (1.83) 2 Weisskopf units (W.u.). 
A good overall agreement is found but for the decay of the state: the B(E2\ 0^" — > 2+) 
value is overpredicted by more than an order of magnitude while the B(E2; 0% — > 2f) value 
is too small by a factor two. This casts doubt on the interpretation of the observed at 
1460 keV as the /3-band head. Previous interpretations of this state are contradictory: it is 
considered as the /9-band head in but as a two-quasiparticle state in |Z| while Casten 
and von Brentano p3| claim it is a collective phonon excitation built on the 7 band. Other 
+ states are found in 166 Er at slightly higher energy p2] but none has the decay pattern in 



agreement with the present calculation. A possible explanation is that collective strength is 
fragmented through mixing with two-quasiparticle states which are absent from the IBM-1 
model space. 

In summary, three-body interactions in the IBM-1 can account for a wide variety of 
7- vibrational anharmonicities in nuclei such as for instance those observed in 166 Er but not 
without substantially changing the moments of inertia of various bands. The knowledge of 
the two double-7 vibrational bands (K n = + and K n = 4 + ) in a single nucleus provides 
a stringent test of nuclear models and, specifically, of the type and strength of three-body 
interactions in IBM-1. More experiments on double-7 vibrations are thus called for since 
they should provide essential information concerning the systematic behavior of these states 
and hence the interactions involved. From the theoretical side, a systematic analysis of all 
three-body interactions and not just those between the d bosons seems in order. Once a fuller 
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knowledge is acquired of the systematic behavior of the interactions necessary to reproduce 
the observed anharmonicities, one may then attempt an understanding on a microscopic 
level. 

We are grateful to F. Iachello, A. Vitturi, and C. Volpe for valuable comments. This 
work has been supported in part by the Spanish DGICYT under contract No. PB98-1111 
and by one IN2P3 (France) -CICYT (Spain) agreement. 
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FIGURES 

FIG. 1. The ratios E? K (as defined in the text) as a function of x- The Hamiltonian (1) is used 
with 9i = 0; the boson number is N = 15. 

FIG. 2. The ratios W K (as defined in the text) as a function of 6i/k for different I. The 
Hamiltonian (1) is used with x = —0.5; the boson number is N = 15. The dashed lines give the 
experimental values for the corresponding ratios in 166 Er. 

FIG. 3. Experimental (a) and calculated (b) spectrum for 166 Er. The theoretical results are 
obtained with the Hamiltonian (1) with k = 23.8 keV, x = —0.55, k' = —1.9 keV, and 64 = 93.9 
keV. The boson number is N = 15. 
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TABLES 

TABLE I. Observed and calculated B(E2) values and ratios for 166 Er. The E2 operator (4) is 
used with e 2 eQ = (1.83) 2 W.u. and x = -0.55. 

B(E2) value or ratio 







Observed 


Calculated 


£(E2:2+ - 


* Ot) (W.u.) 


214 ± 10 a 


214 


£(E2;4+ - 


» 2+) (W.u.) 


311 ± 10 a 


304 


£(E2;2+ - 


» 0+) (W.u.) 


5.5 ±0.4 a 


5.3 


£(E2;0+ - 

S(E2;2| - 
£(E2;0+- 

S(E2;2| - 
B(E2;0+ - 

i?(E2; 2j - 
S(E2;4+ - 

5(E2;2| - 


.2+) 
,0+) 
+ 2+) 
,0+) 

-0+) 


0.49 ±0.19 b 
0.44 ±0.13 b 
3.8 ± 1.3 c (2.2±Jj d ) 
1.3 ±0.4 c (0.9 ±0.3 d ) 


0.21 
6.2 
3.2 
2.5 



a From Ref . @ . 
b From Ref. fFI 



c From Ref. [§]. 
d From Ref. 0. 
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